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ABSTRACT 



Context. Small universe models predicted a cutoff in large-scale power in the cosmic microwave background (CMB). 
This was detected by the Wilkinson Microwave Anisotropy Probe (WMAP). Several studies have since proposed that 
the preferred model of the comoving spatial 3-hypersurface of the Universe may be a Poincare dodecahedral space 
(PDS) rather than a simply connected, flat space. Both models assume an FLRW metric and are close to flat with 
about 30% matter density. 

^ , Aims. We study two predictions of the PDS model, (i) For the correct astronomical positioning of the fundamental 

$H ' domain, the spatial two-point cross-correlation function of temperature fluctuations in the covering space (where 

' the two points in any pair are on different copies of the surface of last scattering (SLS)) should have a similar order of 

I magnitude to the auto-correlation function on a single copy of the SLS. (ii) Consider a "generalised" PDS model 

for an arbitrary "twist" phase <j) £ [0,27r]. The optimal orientation and identified circle radius for a generalised PDS 
model found by maximising relative to in the WMAP maps should yield one of the two twist angles ±36°. 
(N ■ Methods. Comparison of to extends the identified circles method, using a much larger number of data points. 

^ We optimise the ratio of these functions at scales ^4.0/i~^ Gpc using a Markov chain Monte Carlo (MCMC) method 



o 
o 
o 



over orientation {I, b, 6), circle size a, and twist 

Results. Both predictions were satisfied: (i) An optimal generalised PDS solution was found for two differ- 
ent foreground-reduced versions of the WMAP 3-year all-sky map, both with and without the kp2 galac- 
tic contamination mask. This solution yields a strong cross-correlation between points which would be distant 
and only weakly correlated according to the simply connected hypothesis. The face centres are {(Z,6)}i=i_g « 
(3 ■ {(184°, 62°), (305°, 44°), (46°, 49°), (117°, 20°), (176°, -4°), (240°, 13°)} (and their antipodes) to within « 2° ; (ii) This 

QQ ' solution has twist <f) = {+39 ± 2.5)°, in agreement with the PDS model. The chance of this occurring in the simply 

, connected model, assuming a uniform distribution </> £ [0, 27r], is about 6-9%. 

Conclusions. The PDS model now satisfies several different observational constraints. 
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1. Introduction three-dimensional space, e.g. the two-point auto-correlation 

function. 

The past decade and a half have shown considerable growth The predicted cutoff in large-scale power appears to 

in attempts to determine the global shape, i.e. not only the have been confirmed by the first-year observations of the 

curvature, but also the topology, of the spatial comoving Wilkinson Microwave Anisotropy Probe (WMAP) exper- 

section of the Universe, i.e. of the 3-manifold to which a 3- ™ent. With this data, Spergel et al. (2003) pubHshed a 

hypersurface corresponds, informally known as "space" . As ^8^^ approximately equivalent to such a function, i.e. the 

noted by several authors, in particular Starobinsky (1993) black "WMAP data" curve of Fig. 16 in their paper. Their 

and Stevens et al. (1993), a space that is "small" com- figure shows the auto-correlation calculated as a function of 

pared to the surface of last scattering (SLS) cannot contain angular separation, shown against projected spatial separa- 

eigenmodes, which are used for expressing density perturba- tion for a (first year) template-cleaned V map with the kpO 

tions, which are themselves larger than the space itself. This galactic contamination mask. The authors note the surpris- 

should lead to a cutoff of power in statistics representing ingly fiat correlation on large scales and suggest a multiply 

these fiuctuations, above which power should drop to zero, connected universe model to match this function. 

This prediction was made after COBE data were available, In Fig. [Tl we calculated the auto-correlation function 

but before the WMAP satellite was launched. For practical, directly as a function of three-dimensional spatial separa- 

observational reasons, spherical harmonic analyses of tem- tion, not of angular separation, using the 3-year integrated 

perature fluctuations on the 2-sphere are frequently made, linear combination (ILC) map with the kp2 cut. This flg- 

However, a physically more natural statistic to use is one in ure can be approximately compared to the black "WMAP 
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spatial separation in h ^Gpc 



Fig. 1. Estimate of the spatial auto-correlation function, of 
temperature fluctuations in the ILC WMAP map. The estimate 
assumes a simply connected space with Sl^ot ~ 1-0' ~ 0-3' 
using 30,000 points randomly chosen from a uniform distribution 
on the 2-sphere, excluding points falling within the kp2 Galaxy 
mask. Correlations are shown in ^K^ against comoving separa- 
tion in Gpc ranging from zero separation up to the diameter 
of the SLS. An approximate conversion from spatial to angular 
separations on small scales (d <,4/i~^ Gpc) can be obtained by 
setting 19.2ft-^ Gpc = 2 rad = SGOVvr, i.e. Ih'^ Gpc « 6.0°. 



data" curve of Fig. 16 of Spergel et al. (2003), except that 
the present figure shows the auto-correlation calculated as 
a function of three-dimensional spatial separation, not of 
angular separation; the relation between the two is not lin- 
ear. Also, our figure uses the 3-year ILC map with the kp2 
cut, not the 1-year template-cleaned V map with the kpO 
cut. 

In spatial comoving units, we confirm that the auto- 
correlation is very close to flat for separations larger than 
K, 10/i~^ Gpc. The relation between angular and spatial 
scales is, of course, not linear. Equation (fT5|) below (for 
either a spherical covering space or for a flat covering space 
using the limit Rc S> ''SLS — "^Z^) can be used to calculate 
this. 

If the size of the Univers^U is about 10/i~^ Gpc, as 
this figure seems to indicate, then which of the various 3- 
manifolds correctly describes comoving space? Motivated 
by indications that the Universe may have positive cur- 
vature, and using eigenmode-based simulations to study 
the spherical harmonic (C;) spectrum of the WMAP data, 
Luminet et al. (Luminet et al. 2003; Caillerie et al. 2007) ar- 
gue that the Poincare dodecahedral space (PDS) is favoured 
by the WMAP data over an infinite, simply connected 
flat space. Caillerie et al. (2007) state that, by requir- 
ing maximal repression of the quadrupole signal, an op- 
timal total density of f^tot = 1-018 is favoured (for a non- 
relativistic matter density = 0.27 and Hubble constant 
Hq = 70kms~^/Mpc) for the PDS model. Several other au- 
thors (Aurich et al. 2005a, 2005b; Gundermann 2005) have 
also compared simulations for PDS models to the observed 
first-year and three-year Wilkinson Microwave Anisotropy 
Probe (WMAP) maps of the cosmic microwave background 
(CMB). 



^ See e.g. Fig 10 of Luminet & Roukema (1999) for a schematic 
diagram of various definitions of the "size" of the fundamental 
domain, including the injectivity radius, the in-radius, and the 
out-radius. 



In all these studies, both the infinite flat models and 
the PDS models are used in the context of a standard, hot 
big bang model, i.e. where the Universe has a Friedmann- 
Lemaitre-Robertson- Walker (FLRW) metric, perturbed by 
fluctuations that collapse gravitationally to form structures 
such as filaments and clusters of galaxies, and where values 
of the metric parameters consistent with the consensus ob- 
tained during the last decade of observations are adopted: 
the Universe is close to flat on length scales up to the SLS 
and has about 30% non-relativistic matter density. 

In other work, Roukema et al. (2004) used the identified 
circles principle (Cornish ct al. 1996, 1998) to find a specific 
optimal orientation of the PDS model based on the WMAP 
first-year ILC map, and published a tentative set of coordi- 
nates. Key et al. (2007) confirmed the presence of a signal at 
the celestial coordinates, circle radius, and —36° twist pub- 
lished in Roukema et al. (2004), but argue that it should 
be considered a false positive. Using software independent 
of that used in Roukema et al. (2004), updating to the 
3-year WMAP data, and using Gaussian simulations, Lew 
& Roukema (2008) find similar conclusions. They find that 
a local maximum in the statistic used for finding matched 
circles exists for a circle radius ~ 11° and a —36° twist, but 
it is not statistically significant. 

Aurich et al. (2006) also made a circles analysis of 
the WMAP first-year data, using their own estimator and 
weight function, and, in contrast with Roukema et al. 
(2004), Key et al. (2007), and Lew & Roukema (2008), 
did not find any signal at 11°. On the other hand, they 
did find a tentative PDS signal at fitot ~ 1-015, or equiv- 
alently, a circle radius a ~ 40°. The signal is weaker than 
they expected, but the authors note that uncertainties due 
to foregrounds and noise structures in the data make it 
premature to draw firm conclusions. 

A disadvantage of the identified circles approach is that 
it is based on the information in a relatively small number 
of points on the sky map of temperature fluctuations, mak- 
ing it sensitive to small errors in the data or analysis and 
requiring prohibitively long computations. Is it possible to 
generalise from the identified circles principle? 

Moreover, leaving aside the debate about matched cir- 
cles statistics, observably multiply connected models could 
reasonably be said to have satisfied only one prediction so 
far, that of a cutoff in the density fluctuation spectrum on 
a large scale. Can predictions of the PDS model itself be 
tested, for example, using the identified circles principle or 
an extension of it? 

In searches for matched circles, the statistic used is usu- 
ally some variation on what can be considered to be the 
value of the two-point cross-correlation function of observed 
temperature fluctuations at pairs of points in the covering 
space, C(](?'), where the two points lie on different copies of 
the SLS, and the separation of the two points in the pair 
(on the different copies of the SLS) is zero in the covering 
space. 

We can write this function as if it were possible to sam- 
ple temperature fluctuations as point objects at arbitrary 
spatial points throughout the covering space, i.e. including 
arbitrarily low redshifts, as well as points beyond the SLS: 

ec(r) ^ {6T{^^.W[g,{x,,)]).^^^.^^^^ (1) 

averaging over triples (ii , Z2 , j) satisfying d {xi^ , [gj (xi^ )] ) — 
r and gj ^ I (the identity / is removed because we only 
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want the cross-correlation). Here, Xi^ for k — 1,2 are arbi- 
trary points on the SLS considered to be located at their 
comoving spatial positions, gj for j = 1, . . . , 12 are the 12 
holonomy transformations in the binary icosahcdral group 
r = /* operating on the covering space that match 
opposite faces of the fundamental polyhedron of the PDS, 
d{x, y) is the comoving spatial geodesic distance between 
two points x,y in the comoving covering space (i.e. an 
arc-length on the covering space embedded in R^, e.g. 
Roukema 2001). The temperature fluctuations on one copy 
of the SLS, ST{x), are extended to the whole covering space 
via the holonomy transformations, i.e., 

ygjET, dT{x) = 6T[g,{x)]. (2) 

The latter assumption, i.e. Eq. enables rewriting 
Eq. ([1]) in a way that becomes observationally realistic, 

CcM ^ (<5T(a:.J<5T(x,,)),^ ,^^^., (3) 

again subject to the conditions that the average is taken 
over pairs satisfying d{xi^, [5^(2^22)]) = ^'^'^ 9j I- The 
auto-correlation function can then be written 

^p^{r)^{5T{x,,)5T{x.,,))^^^^^, (4) 

averaging over pairs satisfying d{xi^,Xi^) — r. In words, 
^p^ is the 3-spatial auto-correlation function for pairs of 
points on a single copy of the SLS, while is the 3-spatial 
cross-correlation function for pairs of points where the two 
members lie on different copies of the SLS in the covering 
space. 

Now rewrite the statistic S used in matched circles 
searches as follows [see e.g. Eq. (9) of Roukema et al. (2004), 
with the normalisations (denominator and monopole T) ig- 
nored for simplicity], 

S ^ {ST.ST,)^^^^^ 

= ec(0) (5) 

for the case where pairs of points (?, j) are (hypothetically) 
multiply- imaged locations located on a pair of matching 
circles in the covering space. We use r — here because 
that is the defining characteristic of matched circles — a 
pair of matching points is a match because the two points 
are the same space-time points, i.e. they are separated by 
r = when their topological images located on adjacent 
copies of the SLS in the covering space are considered. 

Here, generalise from the zero separation cross- 
correlation £,(j{0) to cross-correlations at larger separations 
r > 0. We can expect that 

?c('^>0)<Cc(0)- (6) 

i.e. that correlations weaken with separation and that, in 
general, ^Q{r) should be approximately equal to the auto- 
correlation function i^^^ ('') on scales much smaller than the 
"size" of the fundamental domain, e.g. the in-diameter 2r_ 
(e.g. fig 10, Luminet & Roukema 1999) 

^Cir > 0) - Uir > 0) , ifr«2r_, (7) 

since there is no statistical, physical distinction between a 
pair of points on different copies of the SLS and a pair of 
points on a single copy of the SLS, apart from effects that 
are not locally isotropic, such as the Doppler effect. 



Since is generally large at small r and small at large 
r, obtaining a large correlation requires using a range of 
length scales r that are relatively small, e.g., 

Cc(r <2/i-i Gpc) ^ ^Air^2h-^ Gpc) (8) 

should be a relatively high positive value if the multiply 
connected model being studied is the correct model. 

Another way of saying this is that this test compares 
the spatial two-point cross-correlation function of mapped 
(in the sense of the holonomy transformation g, which maps 
one copy of a point in the covering space to one of its im- 
ages) and unmapped temperature fluctuations in the cover- 
ing space, where the two points in any pair are on different 
copies of the SLS, against the auto-correlation function on 
a single copy of the SLS. 

To see yet another way of thinking about this, suppose 
that the PDS model at a given orientation and implied cir- 
cle size is correct and that temperature fluctuations occur 
as point objects, each emitting isotropically. The covering 
space cross-correlation function, based on pairs of points 
that are close to one another due to the application of one 
or more holonomy transformations g, but are not close to 
one another on a single copy of the SLS, should then be 
statistically equivalent (apart from the Doppler effect and 
foreground/projection effects) to a sampling from the spa- 
tial auto-correlation function which can be estimated 
using points that lie on a single copy of the SLS. 

This gives us our first prediction to test the PDS hy- 
pothesis: Eq. (0j should only be expected to hold if the phys- 
ically correct PDS (or other) 3-manifold is assumed and is 
modelled at its correct astronomical orientation. If the PDS 
model is correct but a wrong orientation is chosen, then the 
cross-correlation inferred from it on small length scales 
— in the covering space since it is the cross-correlation — 
should represent correlations of pairs of points that in real- 
ity are not close together, due to the erroneous orientation 
that causes mappings from one copy of the SLS to another 
to be incorrect. Since the correlation of distant points is, in 
general, small, we have for this incorrect PDS model: 

Cc(?-;S 2/1-1 Gpc) _ ^^(r-> 2/1-1 Qp^,) 

< eA(^;S 2/1-1 Gpc) (g) 

i.e. 

Cc(?-;S 2/1-1 Qp^) ^ £^p^{r<2h-^ Gpc), (10) 

in contrast to Eq. ((S]). Similarly, if the PDS model is in- 
correct, then arbitrary orientations should also yield small 
cross-correlations on small length scales, as in Eq. pH)) . 

In this paper, the initial aim is to maximise 
^Q(r <2/i~i Gpc) relative to ^^(r<2/i-i Gpc), varying 
PDS models over the parameter space of different orien- 
tations and circle sizes. (In the actual calculations, we used 
a range of scales below and a little above ^ 2h~^ Gpc: see 
Sect. 13.31 and Eq. for details.) This should lead to an 
estimate of the best orientation and circle size for the PDS 
model, given the observational data. 

However, to further test the PDS model, con- 
sider a "generalisation" from the mathematically cor- 
rect PDS model to a class of pseudo-models for which 
the cross-correlation function can be calculated, but for 
which mosl[3 of the members of the class are physi- 

^ The mathematical term "almost all" could be used here, but 
when taking observational and numerical calculation uncertain- 
ties into account, the word "most" is more realistic. 
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cally invalid. Maximising ^Q(r <, 2/i~^ Gpc) (relative to 
^^(r <,2/i~^ Gpc)) over the extended parameter space de- 
fined by this class should then have little chance of yielding 
an optimal model that is valid, unless the PDS model is 
astronomically correct. 

The "generalisation" that we define here is to allow the 
"twist" angle to be arbitrary in [0, 27r]. The twist angle 
can be described as follows. In a single-action, spherical 3- 
manifold (e.g.. Sect. 4.1, Gausmann et al. 2001) thought 
of as embedded in 4-dimensional Euclidean space, R^, any 
holonomy transformation is a Clifford translation that ro- 
tates in about the centre of the hypersphere in one 2- 
plane by an angle of 7r/5, and also by tt/ 5 in an orthogonal 
2-plane. 

From a close-up perspective of the SLS rather than look- 
ing at the whole hypersphere S^, or in other words, from 
a projection into R'^, one 4-rotation can be thought of as 
a translation in M.^ from one circle in a matched circle pair 
to the other member of the pair on the opposite side of 
the SLS, and the second 4-rotation is then thought of as 
a twist around a vector in the translation direction joining 
the centres of the two identified circles. This full motion 
is frequently termed a "screw motion" . The "twist" is the 
second of these two rotation angles. 

Physically, the only two possible twist angles are ±7r/5. 
If the cross-correlation were to be calculated by filling one 
copy of the fundamental domain with a uniform distribu- 
tion of points and then mapping this set of points to copies 
of the fundamental domain in the covering space, then for 
an invalid twist angle, some regions of space near the first 
copy would have either zero, two, three, or more times the 
density of points in the first copy of the fundamental do- 
main. In other words, sharp discontinuities would occur in 
calculating the cross-correlation. 

However, a practical way of estimating the cross- 
correlation is the method represented in Eq. ([5]), i.e. multi- 
plying temperature fluctuations and then averaging them. 
Using this method rather than multiplying densities of 
points, the effect mentioned in the previous paragraph 
should only lead to slightly more or less frequent sampling 
in some regions of the cross product of the covering space 
with itself, not to a modification of the correlations them- 
selves. 

The continuous nature of this method, as opposed to 
the discrete nature of starting with a uniform distribution 
of points in the fundamental domain, can be seen as follows. 
For a given pair of observed locations (xj, Xj) on one copy 
of the SLS, their comoving spatial separation d[xi, g{xj, (j))] 
for a given "generalised" isometry g, which depends on 
the generalised twist </>, changes continuously through the 
values (j) = i^'/Sj not discretely. Moreover, the product 
of the observed temperature fluctuations remains constant. 
Intuitively, we could say that as (j> varies, the geodesies 
formed by pairs [xi, gi{xj , (f))] do not "know" when they 
meet other pairs [xi,g2{xj,(j>)], where gi and g2 are two 
holonomy transformations mapping the fundamental do- 
main to two copies of the fundamental domain adjacent to 
one another in the covering space. 

If the PDS model is correct, then calculating 
£,c{r) using temperature fluctuations and applying 
isometrics using an arbitrary twist angle should give 
^C{r<2h-^ Gpc) < ^j^{r<2h-^ Gpc) as in Eq. 1^, but 
give ^(j{r<2h-^ Gpc) ~ ^p^{r <2h-^ Gpc) as in Eq. ^ 
when the twist angle is correct. 



This gives the second prediction of the PDS model. 
The maximal cross-correlation for this generalised PDS 
model estimated using correlations of temperature fluctu- 
ations should not only exist as a robust maximum, but it 
should also give a twist angle of either ±7r/5. 

If the simply connected, perfectly flat model is cor- 
rect, then the PDS model is incorrect, and it should 
either be difficult to find a robust maximal correlation 
^(j(r ^2h~^ Gpc) in this extended 5-parameter space, or 
else an arbitrary twist angle (j) should result, with only a 
small chance of cj) being close to either of the two values 
expected for the PDS. 

An estimate of how close an "arbitrary" twist angle 
should lie to one of the two PDS values can be made as 
follows. Given the assumption that 

eA(r-> 2/1-1 Qp^) ^ ^^(^j.<2h-^ Gpc), (11) 

we can expect that this arbitrary angle should be selected 
from a uniform probability density distribution on [0, 27r]. 
In principle, depending on the assumptions made about 
the complex statistical properties of the WMAP tempera- 
ture fluctuation maps, it could be possible for this distribu- 
tion to be non-uniform, even for a non-PDS space model. 
However, the estimated spatial auto-correlation function i^^^ 
for r >,10/i~i Gpc is close to zero, as shown in fig 16 of 
Spergel et al. (2003) as a function of angular separation on 
the SLS (S^), and explicitly as a function of spatial sepa- 
ration in Fig. [1] here, so the assumption does appear to be 
supported by the empirical evidence. 

For the uniform distribution assumption, let us define 









, 9tt 




( 




1 




) 



(12) 



for (j) £ [0, 27r]. The chance of the observational optimal 
phase 0\iVMAP being close to 7r/5 or 97r/5 is then 

f 2^ if AS < - 
P (A0 < A0WMAP) = 1 ; M If ^ < < (13) 

The piecewise nature of this function is because; e.g. there 
are four ways in which ip can differ from ±7r/5 by a small 
angle such as 10°, but only two ways in which it can differ 
from ±7r/5 by a large angle such as 100°. 

To search for an optimal solution in the parameter 
space, a Metropolis-Hastings version of a Markov-chain 
Monte Carlo (MCMC) method is used (e.g.. Sect. 4.2, Neal 
1993, and references therein) over the five-dimensional pa- 
rameter space 

{{l,b,e,a,cj,)}, (14) 

where the parameters represent dodecahedron orientation 
(l, b, 9), circle size a and twist phase (f). All parameters are 
initialised as arbitrary angles, except that the circle size is 
constrained to 5° < a < 60° both initially and during the 
whole chain. See Sect. 13.61 for more details. 

We also ran some MCMC chains starting with the pa- 
rameters of the Roukema et al. (2004) "hint" of a PDS so- 
lution with circles of size ^ 11°. If this solution is correct, 
then we would expect the chains to remain localised around 
that same solution. If the solution is wrong and if a correct 
solution exists, then the chains should move towards that 
correct solution, even though starting at a wrong point. If 
the solution is wrong and if no correct solution exists, then 
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we would expect the chains to move randomly and fail to 
find a strong maximum. 

The correlation function definitions are described 
in Sect. 13. 2[ the probability estimator used to com- 
pare multiple-SLS cross-correlation and single-SLS auto- 
correlation functions is presented in Sect. 13.51 and the 
MCMC method is described in Sect. 13.61 Results are pre- 
sented in Sect.m Discussion follows in Sect. [5] and conclu- 
sions are made in Sect. [6l 

For references on cosmic topology in general, please 
see the first known article on the subject (Schwarzschild 
1900, 1998), a short beginner's review of cosmic topology 
(Roukema 2000), more in-depth reviews (Lachieze-Rey & 
Luminet 1995; Luminet 1998; Starkman 1998; Luminet & 
Roukema 1999), workshop proceedings (Starkman 1998; 
Blanlceil & Roukema 2000), and lists of two-dimensional 
and three-dimensional methods (Table 2, Luminet & 
Roukema 1999; Uzan et al. 1999; Roukema 2002; Rebougas 
& Gomero 2004). For background on spherical, multi- 
ply connected spaces, see Weeks (2001), Gausmann et al. 
(2001), Lehoucq et al. (2002), Riazuelo et al. (2004), and 
Luminet et al. (2003). For the identified circles principle, of 
which the present method can be thought of as an exten- 
sion, see Cornish et al. (1996, 1998). 

One correction to a point made in much of the above lit- 
erature concerns the independence of the metric and global 
topology. For example, Luminet & Roukema (1999) wrote, 
"However, general relativity deals only with local geometri- 
cal properties of the universe, such as its curvature, not with 
its global characteristics, namely its topology". The error 
here is that global topology can, at least in certain cases, 
generate a local effect, and thus affect local geometrical 
properties. See the T^, Newtonian weak field limit, heuris- 
tic calculation in Roukema et al. (2007): if the Universe 
contains inhoniogeneities and is not expanding perfectly 
isotropically, then at least in the case, additional accel- 
erations and decelerations exist in the different fundamen- 
tal directions, in such a way that they tend to equalise the 
slightly different expansion rates. 

Comoving coordinates are used when discussing dis- 
tances (i.e. "proper distances", Weinberg 1972, equivalent 
to "conformal time" if c = 1). We write the Hubble constant 
as Ho = 100/ikms-VMpc. 

2. Observations 

The analysis presented here uses the Internal Linear 
Combination (ILCjl all-sky map of the three-year WMAP 
data (Spergel et al. 2007) and the Tegmark et al. (2003) 
foreground-cleaned, Wiener-filtered (T0Hj3 version of the 
same data. (Unless otherwise stated, we refer to the three- 
year data, not the one-year data.) 

For our main analyses, we used either the "kp2" mask 
to eliminate the Galactic Plane (CP) and associated re- 
gions from the analysi^ or no mask at all. We also discuss 

^ | http : //lambda.gsfc .nasa. gov/data/map/ I 
|dr2 /df p/wmap_ilc_3yr_v2 .fits 

http : // space . mit . edu/home /tegmark/| 
[wmap/wiener3yr_map .f its 

^ Data file: http : //lambda, gsfc .nasa. gov7data7map/ 
|dr2/ancillary/w map_kp2_r9_mask_3yr_v2 .f its , map pro- 
jection: ^h ttp : //lambda, gsfc . nas a.gov/product /map/] 

I current / map_images/ f 02_int _mask_b . png j 
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Fig. 2. Schematic diagram of temperature aufo-correlations 
in a simply connected universe, showing "typical" temperature 
correlations STiSTj for two pairs of points (i,j) = (1,2), (3, 4) 
that are widely separated on a single copy of the surface of last 
scattering (SLS). Since the correlation is in general high at low 
separations and low at high separations, the product STiSTj is 
low in both cases. 
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Fig. 3. Schematic diagram of temperature cross-correlations 
in a multiply connected universe (e.g. a Poincare dodecahedral 
space, PDS, model), showing "typical" temperature correlations 
STi STj for two pairs of points (i,j) = (1,2), (3, 4) that are 
widely separated on a single copy of the SLS, as in Fig. [5] The 
holonomy transformation g maps points in space to copies of 
themselves in the covering space. The points xi and g{xi) are 
the same physical point, and similarly, 2:3 and gix^) are physi- 
cally identical. Since the comoving separations d[g{xi),X2] and 
d[g{x3),X4] are both small, and the correlation is, in general, 
high at low separations, the product STiSTj is high in both 
cases. This should only occur if the holonomy transformation 
g is physically correct. This is described statistically in Eq. (|8]). 



analyses with the "kpO" mask. The kp2 mask covers about 
15% of the sky, and the kpO mask covers about 25% of the 
sky. 

We did not carry out any explicit smoothing of these 
maps, though the MCMC method itself could be thought 
of as a method that (when analysed statistically) implicitly 
smooths the information in the data. 
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3. Method 

A GPL (GNU General Public Licence) program CIRCLE^ 
is available for reproducing the analysis as described be- 
low. A typical 12,000 step MCMC chain as described below 
should take about 3-4 days to run on an x86 type processor 
running at about 1.5-2 GHz. 

3.1. Spatial versus angular correlation functions 

Since the naive Sachs- Wolfe effect is due to locally isotropic 
emission, which is what makes the identified circles test use- 
ful for the purposes of testing cosmic topology hypotheses, 
we will assume the same approximation as that required 
for the circles test to reveal a signal, i.e. that the inte- 
grated Sachs- Wolfe effect and the Doppler contribution are 
relatively small on the observational angular scales corre- 
sponding to the length scales of interest. At worst, this 
assumption should weaken any signal found; it is difficult 
to see how it could create a false signal. This assumption is 
written below in Eq. 

3.2. Auto- and cross-correlation functions: vs 

Due to observational constraints, the spatial auto- 
correlation function that we are interested in for this test is 
usually estimated indirectly, e.g. in angle, and calculated in 
terms of spherical harmonic components of the temperature 
fluctuations considered as a function on the sphere. 

However, the angular correlation function on the SLS 
is sometimes calculated. Figure 16 of Spergel et al. (2003) 
shows the main feature useful for this test: to first order, the 
correlation is high at low separations and low at high sepa- 
rations. As mentioned above, the figure also shows both the 
disagreement between the simply connected flat model and 
the observations and the surprising flatness of the correla- 
tion function over most angular scales greater than about 
40-50°. 

Figures [2] and [3] show what is described algebraically 
in Eqs. pO)l and ([8]): in a simply connected model, 
^Q{r<2h-^ Gpc) < ^a(^~2/i"1 Gpc), while in a mul- 
tiply connected model at the correct orientation and iden- 
tified circle size, ^c{r<2h-^ Gpc) - ^js,{r<2h-'^ Gpc). 
For simplicity, these two diagrams show only one holonomy 
transformation gj . A holonomy transformation is an isome- 
try of the covering space used to define the PDS, i.e. gj G F 
where F = /* , the binary icosahedral group, is the group of 
isometrics dividing the covering space into 120 multiple 
copies of the fundamental domain. 

The auto-correlation, is calculated using Eq. (jH). To 
calculate for a given orientation and circle size, we used 
Eq. ([3]), with all 12 of the holonomy transformations (those 
giving the first "layer" of spherical dodecahedra neighbour- 
ing the "first" copy of the fundamental domain) . 

In principle, Eq. ^ could be evaluated over a "large" 
subset of all members of the subgroup generated by the 

^ Version CIRCLES-0.2.4 is the one used for the main 
calculations in this paper, along with CIRCLES-0.2.4.f13.3, 
which differs from the former only in the formatted print- 
ing of galactic longitude to standard output and to a 
file. Various versions of the packet are downloadable from 
http://adjani.astro.uink.pl/GPLdownload/dodec/, The soft- 
ware presently requires medium to advanced GNU/LlNUX, 
Fortran77 and C experience for a scientific user. 




(Surface of Last Scattering) 

Fig. 4. Surface of last scattering and its interior, i.e. the co- 
moving 3- volume of the observable Universe as a subset of S^, 
embedded in for convenience. All lines appearing straight in 
this figure correspond to spatial geodesies in S^, i.e. arcs in R*. 
The triangle to the left illustrates the geometry of a matched cir- 
cle. In the R* embedding, the PDS mapping from one matched 
circle to another rotates by 2n/10 = 36° between circle centres, 
so an arclength of half this (from the observer to one matched 
circle centre), i.e. n/10 multiplied by the 3-sphere radius, Rc, is 
shown. The circle radius as measured on the SLS is shown as a. 
The relevant equation describing this triangle is Eq. (I15|l . The 
right-hand part of the figure shows half the geodesic distance d/2 
between two points Xi,Xj lying on the SLS and the angle Od/2 
subtending this. The relevant equation describing this triangle 
is Eq. (I16|) . The relations between angles and arclengths of the 
sides of these triangles are determined by spherical trigonome- 
try. 



holonomy transformations, e.g. <?25i, 5455512, ■ • ■• However, 
this would be highly redundant. The minimum number of 
calculations that should be done to test a reasonably large 
sample of different pairs would be those within one copy of 
the fundamental domain, plus pairs crossing the boundaries 
of that copy of the fundamental domain. 

This leads to an alternative way of visualising ^q. It 
could be thought of as the auto-correlation function of the 
full set of all temperature "fluctuations" (or pixels for a 
given pixelisation of the sphere) mapped into one copy of 
the fundamental domain, plus those crossing "boundaries" . 

In practice, we find it simpler to work in the covering 
space and to use just the set of holonomy transformations 
that give the neighbours of a single copy of the fundamental 
domain, rather than the full group of 120 holonomy trans- 
formations. 

3.3. Comoving lengths and angular scales 

Physically, the separations of interest are lengths of spa- 
tial geodesies in the covering space — the hypersphere. 
This is used for both the simply connected and PDS as- 
sumptions, since there is no problem having just a single 
SLS embedded in an S'^ covering space. These separations 
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can be thought of as arclengths or the angles that subtend 
them on that hypersphere, of radius Rc (e.g. see eq. (2) of 
Roukema et al. 2004), embedded in Euclidean 4-space (e.g. 
Roukema 2001), but these angles are different from angles 
on the SLS (a 2-sphere of radius t'sls) ^ seen from an 
observer at the centre of that 2-sphere. 

Fig. [5] shows some key relations within the SLS. From 
the left-hand triangle, use of the spherical sine and cosine 
formulae leads to 



I.e. 



tan 



^SLS ^ tan(7r/10) 
Rc cos a ' 



(15) 



for the ratio of the SLS (5^) radius, ''sls, to the radius of 
the full hypersphere {S^)^ Rc, where a is the radius of a 
matched circle (S-^) on the SLS. For the right-hand triangle 
a single application of the spherical sine formula is sufficient 
to obtain an expression for relating the geodesic (3-sphere) 
distance d{xi,Xj) between a pair of points {xi,Xj) to an 
angle 9d{xi,Xj) more familiar to observers of the cosmic 
microwave background: 



. d/2 
sm — — 

Rc 



sm ~^ sm — 
Rc 2 



(16) 



However, what is of interest in the present work is to es- 
timate the spatial correlation functions in comoving space 
rather than angular correlation functions on the SLS. While 
it is possible to define a meaningful relation between d and 
9d when both Xi and Xj lie on the SLS, a more general dis- 
tance d betweeen Xi and gk{xj) for one of the 12 holonomy 
transformations gk cannot easily be interpreted as an angle 
9d, since, in general, gk{xj) does not lie on the SLS. 

For the auto-correlation function, there is only one copy 
of the SLS, so the maximum separation d at which could, 
in principle, be calculated is the diameter of the SLS, i.e. 

d < 2rsLS- 

For the cross-correlation, since copies of the SLS fill (re- 
dundantly, with overlaps) the covering space (the hyper- 
sphere), the maximum separation d that could, in princi- 
ple, be used for estimating is Tr/JcQThis is several times 
larger than 2rgLg. However, correlations across more than 
half the size (e.g. injectivity radius) of the fundamental do- 
main contain information redundant with those at smaller 
separations. 

Moreover, as mentioned above, we focus on a range 
of length scales that are small enough that the auto- 
correlation is (certainly) high and the cross-correlation (in 
the case of a correctly oriented PDS, if it is physically 
correct) is also high. Above we suggested d<2h^^ Gpc. 
Figure [1] shows the auto-correlation function estimated 
from the ILC WMAP map using the kp2 galactic contam- 
ination mask. The correlations are high for d<2h^^ Gpc 
and remain positive for d^Ah~^ Gpc. They are approxi- 
mately zero for separations greater than one SLS radius: 
9.5/1-1 Gpc<d<19h-^ Gpc. 

Since we want to use a range of small scales, not just 
a single scale, the specific range used in calculations here 
was 



5/90 < d/r-gLS < 40/90, 



(17) 



If the covering space were infinite, then there would be no 
upper limit to the distances at which the cross-correlation could 
be calculated, independently of whether or not at these 
separations were meaningful. 



0.5h-^ Gpc <d< 4.0/1-1 Gpc 



(18) 



for iltot ~ 1-0, ^^m ~ 0.3. The fraction 1/90 is by anal- 
ogy with degrees, but as shown in Fig. [4] and Eq. (fT6|) . 
this does not relate linearly to an angle on the SLS and 
is better considered as an arbitrary unit. The approximate 
correspondence in angles on the SLS, as given in Eq. I|16p . 
is approximately 3° ^6d^ 25°. 

Since a relatively small number of points, A'p, is used 
in estimating the correlation, and since we bin into sep- 
aration intervals, the inclusion of pairs with smaller sep- 
arations, d^0.5h~^ Gpc, would be likely to provide only 
a small, noisy contribution to the total correlation. Using 
these closer pairs without introducing too much noise would 
require increasing A^p, thereby slowing down the calcula- 
tions. Nevertheless, this could be examined in future work. 
See Sect. 13.51 for the values of A'p and adopted number of 
bins. 



3.4. Calculating the holonomy transformation from a pair of 
matched circles 

In Roukema et al. (2004), circle pairs were examined 
without calculating holonomy transformations of the 3- 
manifold, since that was not needed. 

Although it is possible to explore the parameter space of 
different possible 3-manifold sizes and orientations directly 
from sets of holonomy transformations without any initial 
reference to identified circles, since the latter are derived 
from the former, we find it easier to start from the latter. 

We represent the parameter space of these different pos- 
sibilities by the coordinates of the face centres of a funda- 
mental domain (spherical dodecahedron), as stated above. 
In this case, it is possible to reconstruct the holonomy trans- 
formation for a circle pair for the PDS corresponding to 
that circle pair, given the positions of the two matching 
circles. Consider the centres of the two corresponding cir- 
cles as euclidean vectors s, Ms in K^, where M is a four- 
dimensional matrix representation of the holonomy trans- 
formation as a left multiplier whose values numerical values 
are not yet known. Using the term "eigenplane" as intro- 
duced in Roukema (2005), these two vectors determine one 
of the two eigenplanes of the isometry M, as described in 
Eq. (15) of Gausmann et al. (2001). 

As in Section 2.3.3 of Roukema (2005), using s and Ms 
in Eq. (19) of Roukema (2005) gives a second vector, t, in 
this same eigenplane, orthonormal to s: 



Ms - (s • Ms) s 
^1 - (s • A-fs)2 ' 



(19) 



where • is the inner product on M*. 

Since Aft lies in the same eigenplane and resolves to the 
same components in the plane as Ms, apart from a rotation 
by 7r/2, it follows that 



Mt = -(t • Ms) s + (s • Ms) t. 



(20) 



We find a third vector u starting with one of the four 
Cartesian axis unit vectors uo (that which is furthest in 
angle from s, t) 



u = Uo - (uo • s)s - (uo • t)t. 



(21) 
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A fourth vector v is defined as the cross product of s, 
t, and u, which can be calculated by writing these three 
vectors as row vectors s', t', and u', respectively and es- 
timating the determinant of the appropriate 4x4 matrix 
using the four othonormal basis vectors in a given or- 
thonormal basis: 



s' 
t' 
u' 

ei 62 e3 64. 



{22) 



Since the zero point of rotation within the u, v plane is 
arbitrary, we then define the vectors Mu and Afv: 



Mu = (cos 0)u + (sin 0)v 
Mv = — (sin 0)u + (cos 0)v 



(23) 



where <j) is the generalised twist phase mentioned in Sect. [T] 
We then calculate M from 



M = [stuv]"i [Ms Mt Mu M^r]. 



3.5. "Probability" estimator 



(24) 



Our primary aim is to see if we can reject the simply con- 
nected model by looking for the signal expected from a PDS 
model. This does not require a true probability for use in an 
MCMC search of parameter space. Instead, it is sufficient 
to have a function similar enough to a probability func- 
tion that the chains will explore the full parameter space 
and converge after a reasonably short amount of computing 
time and that they will show where the maximal correlation 
hes. 

For a given set of parameters {I, b, 9, a, (f>) [see Eq. (fH)) ] 
defining a sixtuplct of circle pairs, A'p = 2000 points are 
selected from a uniform distribution on the 2-sphere. The 
auto- and cross-correlations and are calculated for 
all ^pair = ^^p(^p - l)/2 and A^^p^ir = 12N^ pairs, re- 
spectively, applying Eqs ^ and respectively. Many of 
these pairs, especially in the latter case, fall at separations 
too large to be of interest for our test. We bin the pairs 
into n = 7 bins in the comoving separation range given by 
Eq. (UHl). 

We define the "probability" that the function at a 
given point in parameter space is sampled from the known, 
"true" correlation, i.e. the auto-correlation^^, by assuming 
a Gaussian distribution of errors of width at in each i-th 
bin (unless > ^a(*)) and assuming that the bins are 

independent: 



P(l,b,0,a,> 



n 



[ic(')-eA(')r 




1 + 0.01 



iHci^)<Ui^) 
if ^c{i)>Ui^)- 



(25) 



Values P < 0.01 are replaced by P = 0.01 so that the 
MCMC chains are ergodic (e.g., Sect. 4.2, Neal 1993). 

The definition of the distribution width, cr^ , is based on 
numerical experimentation in order to obtain an optimal fit 
with the MCMC method in a reasonable amount of time: 



(26) 



where A'^; is the number of pairs contributing to the i-th 
bin in the estimate of ^j^. For the largest bin, i = n, which 
gives width fT„ = ^A(n)/2, so a cross-correlation in this 
bin needs to be as high as ^Q(n) > S^p^(n)/2 to contribute 
a probability per bin of exp(— 0.5) « 0.61 or higher to the 

full product. The weighting \J^^ increases (Ji in bins where 

there are fewer pairs (normally those with i < n), thereby 
moderately decreasing the strength of those bins' contri- 
butions to the product in Eq. (|25p . so that bins with rel- 
atively high Poisson errors do not penalise the probability 

too much. In practice, the maximum of ^ reaches ~ 

for A^p = 2000 and n = 7 bins in the range given in Eq. (fT5|) . 

For bins in which ^c(i) > Ca(*)' ^f the upper ex- 
pression in Eq. ([^5]) would lead to penalising these bins for 
giving unexpectedly strong cross-correlations. Would this 
be reasonable? When cross-correlating points that are phys- 
ically close in comoving space but observed at widely sep- 
arated angles, it is conceivable that the avoidance of fore- 
ground effects could lead to in some bins, in 
which case we would want to favour these contributions. 
To know whether any significant effect should be expected 
in practice would require full modelling of the Doppler and 
integrated Sachs- Wolfe effects. 

However, while it is unclear how significant any such 
effect would be, it seems clear that we should not penalise 
the bins in which > "^aC*)- Here, a very slight excess 

probability per bin is applied in these bins, as defined in 
the lower expression in Eq. (|25p . A bin in which the cross- 
correlation is as high as = 2^^(i), which is much 
higher than is likely to occur in practice, would yield a 
probability per bin of 1.01 to the full product, so the con- 
tribution in practice is likely to be small; i.e., there should 
be a slight preference for excess cross-correlations with re- 
spect to auto-correlations. 

3.6. IVIarkov Chain Monte Carlo method: 
Metropolis-Hastings algorithm 

The function P{l,h,0,a,(f) defined in Eq. ([^S]) is used to 
explore the parameter space defined in Eq. p4p . maximis- 
ing P using an MCMC method with a Metropolis-Hastings 
algorithm. 

An informal description of this method follows. Starting 
from a randomly selected point in parameter space, a 
"chain" of successive points in parameter space is calcu- 
lated, such that a new step in the chain is taken with cer- 
tainty when the new probability at a newly chosen test 
step is higher than the present step, but may or may not 
be taken if the probability at the new step is lower. The 
probability of advancing in the latter case is the ratio of 
the new probability to the probability at the present step. 
This enables a chain of steps to seek the region of highest 
probability, while not getting stuck in noisy areas or local 
hills. For a more formal and complete description of the 
method, see, e.g. Neal (1993), and references therein, and 
in particular. Section 4.2 of that paper. 

At each step, the fc-th parameter is chosen randomly 
using a uniform distribution on {1,...,5}. The proposal 
distribution Sk{x,x%) in the Neal (1993) notation is a 
Gaussian distribution on this A;-th parameter, centred on 
the previous state (/, &, 9, a, 0), of width 



^MCMC - 10°, 



(27) 
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where ° indicates great circle degrees. The Gaussian is 
truncated (without renornialisation) in a at a > 5° and 
a < 60°. 

The probability P as defined in Eq. (|25p is used in 
Eq. (4.18) in Neal (1993) for the acceptance function 
A(x,x*). The initial state (point in parameter space) is 
selected randomly from uniform distributions in 

< / < 27r 
< 6 < 7r/3 

O<0< 2tt/5 (28) 
5° < a < 60° 
< < 27r. 

The region of {I, b, 9) space used for choosing the initial 
values is highly redundant. The range in (?, h) defined here 
covers a fraction sin(7r/3)/2 k. 43% of the surface of the 
sphere; while in principle it only needs to cover 1/12 of the 
sphere, i.e., there is a factor of « 5 redundancy. 

The range in the initial choice of "rotation" angle 
is due to the nature of the dodecahedron. A single choice 
of (^, h) defines a class of dodecahedrons that share one 
pair of faces joined by the axis joining (Z, h) to its antipode 
(/ + TT, —b). Rotation about this axis by G [0, 27r/5) gives 
a class of distinct dodecahedra, but rotation by 2i7r/5 for 
any integer i gives the same set of face centres. 

In the chains, the numerical values of (Z, 6, 0) are all al- 
lowed to increase or decrease arbitrarily, beyond the ranges 
chosen initially. This is for calculational simplicity, both for 
speed and for minimising the chance of introducing errors 
into the software. Together, the triple (Z, &, 9) is likely to 
cover only a relatively small fraction of that part of the 
non-periodic Cartesian product space of the three param- 
eters defined by the least and greatest numerical values of 
Z, b and 9 reached in a given chain. 

During analysis, probably the simplest way to reduce 
the redundancy is by converting each triple (Z, 6, 9) to a 12- 
tuplet {(Z, b)i, i — 1, 12}, where the first face centre (Z, 6)i is 
numerically equivalent to (Z, 6), converted to a conventional 
numerical range (e.g. I G [0, 27r), 5 G [—t:/2, 7r/2]). Although 
the set of face centres should, in principle, consist of 6 an- 
tipodal pairs rather than 12 independent sky positions, the 
freedom given to (Z, &, 9) and the stochastic nature of the 
MCMC chains and the sampling used to estimate correla- 
tion functions imply that the two elements of any antipodal 
pair will not in practice be perfectly antipodal. 

The longest dimension in this parameter space is that 
of the phase (p, so this is the dimension that constrains the 
"burn-in" time. The burn-in time is the number of steps 
in the chain needed to find the region of maximal prob- 
ability, i.e. of the strongest cross-correlations relative to 
auto-correlations, independently of where the chain started 
from. Since the longest possible "distance" from the opti- 
mal region (if that exists) is 180°, the number of steps to 
cross this by a random walk can be estimated as 

/ 130° \ 

A^burn-in~ 5 = 1620 (29) 

Vc^MCMC/ 

steps, where the factor of 5 is due to changing only one 
parameter at any step. 

Thus, for our primary estimate of the optimal point in 
parameter space, we use .^chain chains starting with dif- 
ferent random seeds, each run for 12,000 steps. We ignore 
the first 2000 steps and use the last 10,000 steps for further 
analysis. 
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Fig. 5. Full sky map showing the optimal orientation of dodec- 
ahedral face centres based on 100,000 steps in 10 MCMC chains, 
using the ILC map and either the kp2 mask (upper panel) or 
no mask (lower panel), showing face centres for which P > 0.5 
(see Eq. (|25p ). The projection is a Lambert azimuthal equal area 
projection (Lambert 1772) of the full sky, centred on the North 
Galactic Pole (NGP). The 0° meridian is the positive vertical 
axis and galactic longitude increases clockwise. These face cen- 
tres are derived from the MCMC chains without any constraint 
on the twist phase cf>. 

4. Results 

For the main calculation, N^Yiain ~ MCMC chains were 
run, starting with different random seeds for both maps 
(ILC and TOH), with either no mask or the kp2 mask. Each 
run had 12,000 steps, except that for the ILC map and the 
kp2 mask, the total number of chains was greater {N^Yisln ~ 
24). As expected (Sect.[3l), each chain took about 3-4 days 
to run on a ^ 2 GHz x86 type processor on a system running 
a variety of the GNU/Linux operating system. 

The MCMC chains used in this paper 
can be downloaded for independent analy- 
sis from t he file http : //ad j ani . astro . umk . pl| 
/GPLdownload7MCMC/mcmc_RBSG08 . tbz[ 

4.1. Optimal dodecahedron orientation: {l,b,9) space 

FigureslSlIll and[7]show the sky positions {I, 6)i=i,i2 implied 
by the (Z, b, 9) triples in the MCMC chains for which P > 
0.5 (see Eq. ([25)) ). For consistency between the plots, only 
10 of the 24 available chains are used to create the plot for 
the ILC, kp2 case, even though for calculations, we used 
all 24. The information in the SGP projections (Fig. [6|) is 




Fig. 6. Full sky map showing the optimal dodecahedral face 
centres for the ILC map with the kp2 mask (upper panel) and 
no mask (lower panel), as for Fig. [S] but centred on the South 
Galactic Pole. The 0° meridian is the negative vertical axis and 
galactic longitude increases anticlockwise. 



Fig. 7. Full sky map showing the optimal dodecahedral ori- 
entation for the TOH map with the kp2 mask (upper panel) 
and no mask (lower panel) , as for Fig. [S] centred on the North 
Galactic Pole. 



equivalent to that in the NGP projections (Fig. [5]), so we 
do not show further SGP projections. 

It is clear that a consistent, preferred dodecahedron ori- 
entation exists for both the ILC and TOH maps, both with 
and without masking for contamination by residual fore- 
grounds with the kp2 galactic contamination mask. 

We have estimated the preferred values of {I, b)i for the 
ILC map using the kp2 mask, for which we have 24 chains. 
Since any individual MCMC chain may fail to find the 
global maximum due to starting at a distant point in pa- 
rameter space and random walking in the wrong directions, 
if we search for preferred values of (l,b)i separately in each 
chain and then average them and find the dispersion in the 
estimates, we risk having a large dispersion due to individ- 
ual chains. On the other hand, if we search for preferred 
values of {I, b)i in a concatenation of all 24 chains together, 
we will not have a straightforward way of estimating the 
uncertainty in the estimates of {l,b)i. Hence, we compro- 
mise by grouping the 24 chains together into four groups of 
six chain^. 

For a given group, steps 2001 to 12,000 from each of 
the six chains are concatenated. An approximate (to within 
~ 5-10°) estimate of the optimal (Z, b)i values evident in 



FigsliHZlis 

{(/,fe),}o = {(5°, 5°), (50°, 45°), (120°, 25°), 

(180°, 60°), (245°, 20°), (315°, 50°)} (30) 

and their six antipodes^ This initial approximation is iter- 
ated to yield a more precise estimate. Each iteration uses 
the mean values of points lying within an angle (3j of the 
(j — l)-th estimate, starting from j ~ 1, where j = rep- 
resents the intial approximation. The angular radii start at 
/3i = 30° (so that most of the sphere is covered), decrease 
by 1° for the next 10 iterations, and then remain constant 
at /3j>ii = 20°, covering roughly half of the sphere, until 
the iteration for a given face number converges. For a given 
threshold fmin; only points for which P > Pj^in are used. 
The means are calculated as means of Cartesian {x,y,z) 
vectors formed from the individual sky positions {I, b)i. 

These four groups of six chains are then treated as 
four independent estimates in order to estimate the un- 
certainties due to our MCMC estimation method. These 
uncertainties are standard errors in the mean calculated 
in Cartesian space, converted to a one-dimensional uncer- 
tainty in great-circle degrees assuming (conservatively) that 



The ordering is random, defined by the numerical values of 
the random seeds used for the different chains. 



^ For example, these values can be read off a Mercator pro- 
jection of the same data. 
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Table 1. Sky positions of the best estimate of the six face 
centres for the fundamental dodecahedron for the ILC map 
with the kp2 mask, as shown in Fig. [5l for three values of 
the minimum probability P^i^i (^^^ Sect. |4?T|) . 



Table 2. Estimates of matched circle radius a and twist 
phase (f) using the same points in the MCMC chains used for 
the optimal dodecahedral face solution indicated in Table [1] 
(see Sect. |4?2l regarding systematic uncertainties in estimat- 
ing a). 
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Fig. 8. Distribution of a and states in the MCMC chains 
of the dodecahedral solution used in Table [1] A single point is 
shown independently of whether the chain spends a long time 
at that point — implying that it is a highly probable point — 
or a shorter time at that point, so the density of points does not 
fully represent the information in the distribution. In addition 
to the main cluster of points with ~ +36° (see Table [2] for a 
precise estimate), a secondary weak cluster of points exists for 
—36° = 324°, and an even weaker, tertiary cluster exists 
for (h ~ 280°. 



most of the uncertainty approximately lies in a single di- 
mension: 



over iV = 4 groups of chains. 



(31) 
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Table 3. Estimates of matched circle radius a and twist 
phase (j) for the secondary and tertiary weak features visible 
in Fig. 
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Apart from restricting (j> to the ranges 87r/5 < (p < 2-k and 
67r/5 < (j) < Sn/b respectively, the method is the same as 
for Table [g 




Fig. 9. Schematic diagram showing the approximately parallel 
positioning of two "approximately matched annuli" or "approx- 
imately matched discs" in two copies of the SLS, as per Fig. [l] 
when the matched circle angular radius a is relatively small. 
The centres of the two matched circles projected onto their re- 
spective spheres (copies of the SLS), i.e. P and P', are separated 
by 2[rsLS " (^/lO) Re]- 



The resulting numerical estimates are listed in Table [TJ 
This table lists sky positions of the best estimates of the 
six dodecahedral face centres for the ILC map with the 
kp2 mask, as shown in Fig. [S] The columns show minimum 
probability ^'rninj f^*^^ number i, number n of MCMC steps 
contributing to the estimate obtained from the final itera- 
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tion, galactic longitude I and latitude &, and the standard 
error in the mean between these four estimates in great 
circle degrees, cr((;^b)) [see Eq. ([3T|) ]. The other 6 faces are 
directly opposite with identical errors. As in Fig. [5l these 
face centres were derived from the MCMC chains without 
any constraint on the twist phase cj). 

The final iterated solutions listed in this table are very 
insensitive to the initial approximation: varying the latter 
by up to ~ 10° leads to identical final iterations. The table 
shows that differing thresholds P^[j^ give slightly different 
solutions, especially in the coordinates of the 5-th face, but 
these are consistent with each other within the estimated 
uncertainties. 



4.2. Circle size a 

The optimal dodecahedral face solution shown in Figs (SHU 
and Table [T] is inferred from the data without regard to 
whether the MCMC chains have converged to a preferred 
circle size a and/or twist phase (p. 

We estimated the latter using the same points in the 
MCMC chains used for the estimate in Table [TJ For the 
circle size a, since a is restricted to the range 5° < a < 60°, 
we take the mean and standard error in the mean of a as 
a scalar value. 

However, there are some reasons to be prudent about 
the estimate of a obtained here. Firstly, relatively large 
changes in a for a fixed dodecahedral face set {I, b, 9) and 
a fixed twist (j) can lead to relatively small changes in the 
distances between a given pair of "close" points on two dis- 
tinct copies of the SLS, especially when a is "small". This 
is related to the reason why we expect the cross-correlation 
method presented in this paper to work: there should be not 
only many "perfectly matched" pairs of points, in the sense 
that they are separated by a comoving distance of zero on 
a single copy of the fundamental domain, but there should 
also be many other close pairs, extending the matched cir- 
cle to a "matched annulus" , or in the case of a small circle, 
a "matched disc" . 

Figure [9] illustrates this. For a fixed twist (f>, each pair 
of points on the two parts of the sky inside of two matched 
circles of radius a has a separation of at most the separation 
dc{a) between the two centres of the circles P and P', as 
projected onto the copies of the SLS (not the circle centres 
in From Fig. O and Eq. ([15]) we have 



4(a) = 2[rsLs 



2Rc atan 



(vr/10) Rc] 

tan(7r/10)\ tt 
cos a / 10 



(32) 



For r^tot > 1.01, a circle radius as large as a = 35° gives 
dc < 4:h^^ Gpc. Hence, given that we correlate pairs with 
d^A.Oh^^ Gpc, the main cluster of points in Fig. [8] could 
be interpreted as a "true" matched circle size of ~30-40°. 

Secondly, since the width of the step size in the MCMC, 
'''MCMC — 10° [Eq. P7)) ]. is a large fraction of the range 
in a, it is possible that the chains did not have enough 
freedom to favour a value close to either of these limits, so 
the statistical nature of the precision of any estimate in a is 
less well-established than for the other parameters, where 
the chains are totally unrestricted. 

Inspection of Fig. [8] suggests that the upper hmit, 
a < 60°, did not constrain the chains. However, there is 



a sharp cluster of points at a = 5°. This is partly an arte- 
fact due to not renormalising the Gaussian used in deciding 
whether to move to the new step in an MCMC chain at the 
lower cutoff of a = 5°, as mentioned in Sect. 13.61 Lack of 
renormalisation artificially increases the probability (by up 
to a factor of two) of staying near this boundary. This clus- 
tering may also be caused by the increased noise at small 
circle radii. In any case, if the true value of a were close to 
the a = 5° limit, then the distribution of a values around 
this true value as estimated by this MCMC method could 
not be a two-sided Gaussian distribution, and even a one- 
sided Gaussian distribution would not necessarily occur. 

Both of these caveats should be kept in mind when using 
the value a « 21 ± 1° from Table [2| the systematic error in 
this estimate is likely to be much larger than the random 
error. In this table, as in Table [U the standard errors in 
the mean are obtained by treating the four groups of six 
MCMC chains as independent tests. Listed are minimum 
probability i^niiiu number n of MCMC steps contributing 
to the estimate, a and its standard error in the mean a(^a) j 4' 
and its standard error in the mean cr(0), all four in degrees. 
The number n can be a non-integer since for a given MCMC 
step, it is possible that some of the face centres fall within 
the convergence radius of the final iteration as described in 
Sect. 14. 1[ but other face centres do not. 

4.3. Twist phase (p 

As for a, we estimated using the points in the MCMC 
chains used for the optimal dodecahedral face solution pre- 
sented in Table [1] Since (j> is free to decrease or increase 
arbitrarily in the MCMC chains, we convert (j) values for 
these individual points to Cartesian [x, y) coordinates be- 
fore averaging and reconvert after averaging. The standard 
error in the mean is conservative, assuming that all the 
error occurs in one dimension. 

Table [2] shows the estimate of 4> for slightly differing 
thresholds P^am- Clearly, the value is close to -1-36° within 
about one standard error in the mean, with only a slight 
dependence on the choice of Pniin- Figure [8] shows the dis- 
tribution of a and 4> values corresponding to the favoured 
face centre orientation. 

Only one point is displayed for any given state, inde- 
pendently of how long a chain spent at that point, so the 
visual density of points does not fully match the statistical 
distribution. Hence, the mean values listed in Table [2] can 
be offset from the visual centroids in the figure. 

What is clear in the figure is that, in addition to the 
twist phase 4> ~ -1-36°, a weak secondary feature and a very 
weak tertiary feature are present, of which the first of these 
is consistent with the interesting angle of (/) « —36°, to 
within the estimated uncertainty. 

We estimate the parameters of these two weak fea- 
tures by using a subset of solution points arbitrarily cut 
at 87r/5 < cj) < 27r for the secondary feature and 67r/5 < 
(j> < Sir/ 5 for the tertiary feature. The results are listed 
in Table [3] Comparison of the numbers of steps n in the 
MCMC chains contributing to the different signals, as listed 
in Tables [5] and [31 shows that the secondary and tertiary 
signals are represented by about 10 and 30 times less steps 
(respectively) than the primary signal. 

Clearly, the secondary feature is within l-2cr (standard 
error in the mean) of — 7r/5, with approximately the same 
estimated mean circle radius a « 21 — 22° as the primary 



Roukema et al.: Optimal phase of PDS for WMAP sky map 



13 



I.Or 




Fig. 10. Full sky map showing the optimal dodecahedral orien- 
tation for the ILC map with the kp2 mask, as for Fig. [S] centred 
on the North Galactic Pole, from an MCMC chain starting at 
the PDS orientation and circle size suggested in Roukema et al. 
(2004), for an initial twist of (/> = — 7r/5. The optimal orienta- 
tion is clearly very close to what is found from arbitrary initial 
positions, shown in the previous figures. 




5 10 15 

spatial separation in h~^Gpc 

Fig. 11. As for Fig.[T] together with an estimate of the spatial 
cross-correlation function, (thick curve), made assuming the 
Poincare dodecahedral space with the best estimate parameters 
listed in Tables [T] and [51 corrected to an exact dodecahedral 
solution: {l,b,e,a,(l)) = (184.0,62.0,34.0,21.0,36.0) in degrees. 
Values of are meaningful up to 9.0h~^ Gpc, half the dis- 
tance from the centre of one face of the fundamental domain to 
the matching face. 



signal. Only one of the two twist radii can be valid for the 
spatial matching of opposite faces of the fundamental do- 
main of a PDS model, and they cannot both be valid simul- 
taneously. On the other hand, the density fluctuations in a 
PDS model must resolve into the eigenmodes of the PDS, 
not of infinite flat space. It is conceivable that some sort of 
harmonics could occur, so that secondary and tertiary weak 
features exist in addition to the main cross-correlation. 



4.4. Testing the Roukema et al. (2004) hypothesis 

If the PDS solution found here is correct, then the solution 
suggested in Roukema et al. (2004) cannot be correct (in- 
dependently of any other arguments), since the parameters 
are very different. 

To test this explicitly, we ran 10 MCMC chains on the 
ILC map with the kp2 mask, starting at the Roukema et al. 




184.0 62.0; 4.0 -62.0 
5 10 15 

distance along circle (h~^Gpc) 



Fig. 12. Matched circle pair 1: this and the following fig- 
ures show matched circles for the best PDS model found here, 
corrected to an exact PDS solution as indicated in Fig. 1111 
Temperature fiuctuations are from the ILC 3-year WMAP map 
using the kp2 galactic contamination mask, shown in mK against 
comoving distance along a circle in Gpc. The coordinates 
of the two dodecahedral face centres {l,b)i are indicated. Solid 
lines show the northern galactic member of a pair; dashed lines 
show the southern member. The pixels used for these plots con- 
tributed almost nothing to the method used for finding the opti- 
mal cross-correlation, since pairs of points at nearly zero implied 
spatial separation are rare. (See Sect. 14.71 ) 



(2004) solution. If the latter solution were correct (and if 
this method is correct), then these chains should remain 
around that solution. In fact, as Key et al. (2007) appear 
to agree, the Roukema et al. (2004) solution is at least a lo- 
cal maximum for correlations along matched circles (though 
they argue that it is not statistically significant, as do Lew 
& Roukema (2008) using a different method), so the MCMC 
chains could conceivably remain "stuck" at this local max- 
imum even if a global maximum exists elsewhere in the 
parameter space. However, Fig. [10] shows clearly that the 
MCMC chains failed to remain near the Roukema et al. 
(2004) solution and moved instead to the one found in this 
paper. 

4.5. Probability in the case of a simply connected model 

What is the chance that a twist close to -f7r/5 could have 
occurred in a simply connected model? Conservatively tak- 
ing the 2a upper limits of the estimates of the twist phase of 
the global maximum found by the MCMCs listed in Table[21 
i.e. {(p — 7r/5| <, 5.6-7.8°, and using the probability distribu- 
tion for a simply connected model represented by Eqs (jl2p 
and p^ . the chance of finding (j) to be this close to one 
of the two PDS values in the case of a simply connected 
model is about 6-9%. 



4.6. Cross-correlation of the preferred solution 

Figure [TT] shows the spatial cross-correlation function es- 
timated by assuming that the solution found here is correct 
(thick curve), together with the auto-correlation function 
calculated assuming that space is simply connected. The 
plot is for our solution corrected to an exact dodecahedral 
solution as indicated in the figure caption. 

It is clear that is well above zero in the range 
O.bh^^ Gpc <, d <, 4.0/i"^ Gpc and is of about the same or- 
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Fig. 13. Matched circle pair 2, as in Fig. [121 



Fig. 16. Matched circle pair 5, as in Fig. 1121 
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Fig. 14. Matched circle pair 3, as in Fig. \T2\ The kp2 galactic 
contamination mask cuts severely into these two circles. 



Fig. 17. Matched circle pair 6, as in Fig. [T2] 
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Fig. 15. Matched circle pair 4, as in Fig. [12] As in Fig.[T4l the 
Galaxy cut is severe. 



der of magnitude as the auto-correlation in fact, it is 
above zero in the range 6.0h~^ Gpc. This is clearly con- 
sistent with the PDS hypothesis. On the other hand, ac- 
cording to the simply connected hypothesis, these strong 
cross-correlations are correlations of temperature fluctua- 
tions at points located at positions on the sky separated by 
well above 10h~^ Gpc. Given the nearly flat, zero auto- 
correlation measured on scales above 10h~^ Gpc, these 
temperature fluctuations should be very weakly correlated 
according to the simply connected hypothesis. 



What is striking in Fig. [TT] is that finishes — at 
half the distance separating matching faces of the funda- 
mental domain — just where seems to become flat and 
zero-valued. This is qualitatively consistent with a PDS in- 
terpretation of the "lack of power on large scales" . 

4.7. Implied matched circles 

The method used in this paper uses random pairs of points 
of which each point is selected randomly (uniformly) on 
the 2-sphere, focussing especially on pairs of points that 
are relatively close. 

However, since we used a weighting to focus on bins 
in pair separations in a way that favours bins with more 
pairs (Sect. 13. 5[) . and since we used d<4:h~^ Gpc, corre- 
sponding to 6d^25° [Eq. (US])], most of the correlations 
used to detect the optimal solution are from pairs with 
d 3> 0.5/i~^ Gpc, i.e. 6d ^ 3°, well above our short separa- 
tion cutoff and well above the map resolution. 

While our match is optimal on the scale of several 
h^^ Gpc, it does not necessarily imply a good set of 
matched circles using individual temperature fluctuations. 
Moreover, the caveats on the circle size, explained in 
Sect. I4.2i suggest that, if the PDS solution found here is 
correct, then its "perfect" matched circle size may well be 
rather different from the best estimate found here. On the 
other hand, the first of these caveats suggests that circles 
should still be approximately matched for wrong estimates 
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Fig. 18. Full sky map showing the optimal dodecahedral ori- 
entation for the ILC map with the kpO galactic contamination 
mask, as for Fig. O centred on the North Galactic Pole. 
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Fig. 19. Full sky map showing the optimal dodecahedral ori- 
entation for the TOH map with the kpO galactic contamination 
mask, as for Fig. O centred on the North Galactic Pole. 
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Fig. 20. Full sky map showing the optimal dodecahedral ori- 
entation for f'jjjjji > 0.5, as in Fig. [5] centred on the North 
Galactic Pole (upper panel), and the corresponding a, (j) values 
in degrees (lower panel) for the WMAP Q frequency band map, 
in which the Galaxy very strongly dominates, with no galactic 
mask. Clearly, the signal is very different from what is found in 
the foreground corrected maps. 



of a, as long as the true a is "small" (e.g. ^35°) and/or 
the wrong estimate is not too far from the correct estimate. 

In Figs [T2l - [T7l we show the matched circles for the six 
circles for the corrected solution used in Fig. [TTJ The two 
circles in pair 1 clearly have similar overall shapes to one 
another, as do the two circles in pair 6. Circle pair 2 has 
some regions that match and some of which do not match. 
Circle pairs 3, 4, and 5 are badly affected by the kp2 cut, 
but several of the uncut regions do seem to show some 
correspondence. 

The severe cuts due to the kp2 mask in circle pairs 3, 

4, and 5 illustrate a difference between a search for "pure" 
identified circles and the cross-correlation method, as pre- 
sented in this paper. Since the cross-correlation method is 
not limited to zero separation pairs (in fact, it uses very 
few pairs close to zero separation), it has the advantage of 
being able to cross-correlate "moderately" close pairs and 
obtain some signal for the holonomy transformation g cor- 
responding to a severely cut circle pair, despite the severity 
of the cut. 

5. Discussion 

The PDS model has now satisfied several predictions: not 
only is there a large-scale cutoff in power, but now we find 



that (i) a PDS solution maximising the cross-correlation on 
implied small scales relative to the auto-correlation on cer- 
tainly small scales exists, and the implied small-scale cross- 
correlation is of a similar order of magnitude to the auto- 
correlation on certainly small scales, and (ii) the favoured 
twist of the "generalised" PDS model agrees surprisingly 
well with one of the two possible twists to within the accu- 
racy of the method used. Does this mean that the simply 
connected model should be rejected in favour of the PDS? 

Focussing on the second of these tests — the 6-9% 
chance in the case of a simply connected model of having a 
twist as close to ±7r/5 as what we find here — there are nu- 
merous astrophysical, mathematical, and software caveats 
to be considered. Among these are: 

(i) Could this be an effect of galactic contamination that 
happens to mimic a dodecahedral symmetry with a pre- 
ferred +tt/5 twist when matching dodecahedral faces? 

(ii) Could modifications to the infinite flat model such as 
our location in a local void (e.g. Inoue & Silk 2006) be 
sufficient to modify Eq. so that « 7r/5 is more 
likely? 

(iii) Could the nature of the comparison test between the 
auto-correlation and cross-correlation functions itself 
favour a ±7r/5 twist? Could the MCMC algorithm it- 
self be at fault? 
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Table 4. Estimates of matched circle radius a and twist 
phase (j) as per Table [21 but for the ILC map without any 
mask, and the TOH map with and without the kp2 mask, 
each based on only N = 2 "independent experiments", i.e., 
two independent concatenations of five MCMC chains. 



min 




9 






ILC no mask 






0.4 


6103.5 18.7 0.9 


35.1 


0.7 


0.5 


3513.2 17.6 0.5 


37.4 


0.2 


0.6 


1735.0 17.4 0.4 


38.3 


0.6 




TOH kp2 






0.4 


7144.3 22.1 0.8 


30.8 


4.3 


0.5 


4139.3 21.9 0.4 


30.9 


3.6 


0.6 


2173.8 22.9 0.3 


28.7 


5.3 




TOH no mask 






0.4 


7909.5 19.6 0.4 


27.5 


4.1 


0.5 


4549.8 19.3 0.3 


30.8 


2.4 


0.6 


2203.2 19.1 0.2 


31.2 


1.8 



(iv) Does a simply connected, infinite flat model really im- 
ply a uniform distribution for twist phases 0? 

(v) The angle 7r/5 as a rotation angle in M"* from one copy 
of an SLS to another copy of the SLS is necessarily built 
into the software. Although this rotation is orthogonal 
to the twist, could it be possible that, due to a sub- 
tle programming error, this angle is also built in as a 
favoured twist angle? 

Alternatively, numerous other checks can be made to 
test the PDS interpretation of this result, including: 

(vi) Can this solution be confirmed on smaller scales? 

(vii) Are the "secondary" and "tertiary" solutions listed in 
Tables |3] consistent with PDS harmonics (eigenmodes)? 

(viii) Is this solution compatible with other PDS analyses 
of the WMAP data? 

5.1. Galactic contamination? 

Could the dominant signal be an effect of galactic contami- 
nation? Figures [SH7] show that whether the kp2 mask, cov- 
ering about 15% of the sky, or no mask at all is used, the 
same solution is clearly dominant, both for the ILC map 
and the TOH map. If galactic contamination were signifi- 
cant, it should show up more strongly in the plots with no 
masking. Indeed, Figs [5]-[7] show that some weak amount 
of extra signal, offset from the main solution, does seem to 
appear when no mask is used, consistent with an interpreta- 
tion as galactic contamination. The main solution remains 
dominant. 

These figures only show the orientation of the dominant 
solution. Could the twist itself be due to the shape of the 
mask, e.g. due to correlating "western" mask edges with 
"eastern" mask edges? This seems unlikely, but Table |3] 
confirms that the twist is present with very similar values 
when no mask at all is used. Each case in this table is for 
10 chains. These are divided into two groups, each of five 
chains. The estimates for a given map/mask combination 
are thus based on only N = 2 "independent experiments" , 
each of five chains concatenated, rather than iV = 4 groups 
of six chains, so they should be less precise than those for 
the ILC kp2 case, hsted in Tabled 



The estimates in this table are for smaller numbers of 
chains and so less precise than those for the ILC map with 
the kp2 mask, but they are clearly consistent with the lat- 
ter. The signal is therefore unlikely to be an artefact due 
to the shape of the mask. 

Curiously, the best estimates of (f>, together with their 
uncertainties cr^^^ in Table [3] for the TOH map with and 
without the kp2 mask, seem to indicate that the further 
the mean is from 36.0°, the larger the uncertainty cr^^^. We 
have not attempted to interpret this statistically, though it 
does hint that an alternative way of finding best estimates 
could lead to a best estimate of (j) closer to 36.0°, with a 
smaller uncertainty. As stated at the beginning of Sect. 01 
our MCMC chains have been made publicly available for 
independent analysis. 

Returning to the question of possible galactic contam- 
ination, could it be possible to be even more conservative 
and mask more of the sky while not losing the signal? The 
two-point correlation function depends on the numbers of 
pairs of points available for sampling a domain. Reducing 
the coverage of sky by a fraction / implies a reduction in 
the number of pairs by (1 — /)^. 

This can also be thought of in terms of matched or 
nearly matched circles. For most choices of the twist (/), 
the masked points in one member of a circle pair will cor- 
respond to unmasked points in the other pair. Obviously, 
such pairs cannot be used, even though one of the points lies 
outside of the masked region. Hence, the masking of pairs 
grows faster than the masking of individual pixels. The se- 
vere cuts in about half of the circles shown in Figs [T2HT7I 
illustrate this. 

Any further reduction in the amount of sky used risks 
increasing the presence of noise. A stronger mask used for 
some studies is the kpO mask. This covers about 25% of the 
sky; i.e., about 46% of pairs of points used for estimating 
correlations are masked. Figures [TSl and fT9l show the {l,b)i 
estimates implied by the optimal {I, b, 9) values in the ILC 
and TOH maps, respectively, for the kpO mask. Both figures 
contain the solution found with no mask or the kp2 mask, 
but an additional set of optimal points is present, weakly 
for the TOH map and strongly for the ILC map. 

Which of these is more correct? Tegmark et al. (2003) 
used a method expected to contain less non-cosmological 
signal from foregrounds and detector noise outside the CP 
than the ILC method. This would imply that the additional 
signal that is strong in Fig. [TSland weak in Fig. [12] is more 
likely to be noise than signal. 

In their own study of the first-year ILC and TOH maps 
by an identified circles method, Aurich et al. (2006) find in 
their Fig. 9 that that the TOH map gives a sharper signal 

— at r^tot ~ 1.015 — than the ILC map. While this could 
be a coincidence, it could also indicate that the TOH map 
correctly reduces foreground and detector noise contamina- 
tion in such a way that two quite different methodologies 

— that of Aurich et al. (2006) and our own — yield sharper 
signals favouring a PDS model. 

Given that there is disagreement between the two differ- 
ent maps regarding the additional signal, but the main sig- 
nal remains present, galactic contamination does not seem 
to be a strong contender in explaining the main signal found 
for either no mask or the kp2 mask. 
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Fig. 21. Schematic derivation of the preferred twist of a 
matched circles pair implied by a map dominated by the 
Galactic Plane (GP). The SLS is shown by a circle. The GP 
is shown as a horizontal edge-on, solid disc. The sky map can 
be thought of as a binary map, zero everywhere except at the 
GP where it is positively valued. One copy, AB, of a matched 
circle is shown to the left, at an arbitrary angle with respect to 
the GP, in projection, containing point A, far from the GP, and 
B, intersecting with the GP. The circle intersects with the GP 
at only two points, B and another point behind B in this pro- 
jection. Dashed arrows indicate the translation (with no twist) 
from one side of the sky to the other, so that A maps to A', B 
maps to B'. See Sect. 15.01 



5.1.1. The expected galactic signal 

Another test is to see what signal is expected from the 
Galaxy. We ran four 12,000 step MCMC chains on the Q (41 
GHz) band (smoothed) WMAP maif°l. This map is clearly 
dominated by the Galaxy. Figure [20l shows the "preferred" 
maximal cross-correlations for Pjjiin > 0.5 and the corre- 
sponding circle sizes and twists 0. 

Since the Galaxy signal is nearly axisymmetric, the 
preferred face centres (/, 6)i vary approximately uniformly 
across the whole sphere, with a strong axisymmetry around 
the north-south axis, and some small variation with lati- 
tude. This is not at all similar to the dominant signal found 
in the ILC and TOH maps, which at best could be seen as 
approaching some sort of axisymmetry only by shifting one 
of the face centres to the NGP (or SGP). 

The lower panel in Fig. [20] shows a very strong prefer- 
ence for (j) = TT. This does not match the dominant signal in 
the ILC and TOH maps. Moreover, a favoured twist phase 
of (f) ~ TT should be expected from a map dominated by a 
strong signal approximating that of the GP, for (nearly) 
arbitrary face centres. 

This is shown schematically in Fig. [^T] Consider a sky 
map that is zero everywhere except on the GP, where it has 
some positive value. A circle AB that is translated with zero 
twist to the other side of the sky, i.e. to A'B', fails to match 
the circle A'B'. B is close to the GP but B' is not close to 
the GP, and A is far from the GP, while A' is close to the 
GP. Keeping in mind that the sky here is S'^, it is clear 
that a twist of 180° is required so that the two points that 
intersected with the GP in the first matched circle again 
intersect the GP after the translation. In other words, the 
favoured twist angle is 180°. This geometry applies for arbi- 
trary circles, except for those nearly parallel or orthogonal 
to the GP. Hence, the signal shown in the lower panel of 
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Fig. 22. Comparison of simulated "generalised" PDS maps 
with random input twists (j>^, to the estimated (output) twists 
(j)o, both in degrees. The two correlate well, with no sign of any 



tendency to favour ±36°. The r.m.s. difference between 
(l>o is 20.2°. See Sect.[53]for details. 



and 
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[dr2/skymaps/3yr/wmap_band_smth_imap_r9_3yr-Q-v2 . f its| 



Fig. [20] is just what can be expected from a galactic signal, 
and does resemble our detected signal. 

If the signal we found is caused by the Galaxy, then it 
is not due to the dominant part of the Galaxy signal, but 
rather to a more subtle aspect of the galactic signal which 
by chance has a PDS-like symmetry. 

5.2. Modifications to the infinite flat model 

Is there any modification to the simplest version of the in- 
finite flat model, such as modelling our location in a lo- 
cal void (e.g. Inoue & Silk 2006) or other methods of at- 
tempting to explain the lack of large-scale power or the 
"Axis of Evil" (e.g. Land & Magueijo 2007, and references 
therein), that would explain a dodecahedrally symmetric 
pattern of twists by +7r/5? This could be possible, in prin- 
ciple, although in that case, the question of which model 
is preferred by Occam's Razor would arise. This topic goes 
beyond the scope of this paper. 

5.3. Is 4> ~ ±7r/5 built into our method or software? 

Could the results found here be somehow built into the 
geometrical, correlational, or other algorithmical aspects 
of our method, or be caused by one or more bugs in the 
software? Extensive cross-checking of different parts of the 
CIRCLES package indicates no sign of an error significant 
enough to create an artificial signal similar to the one found 
here. However, in principle, any moderately complex piece 
of software can always contain bugs. 

As discussed in Sect. I5.1.1[ analysis of the WMAP Q 
map gave a clear signal consistent with what is expected 
— an axisymmetric, nearly uniform distribution of pre- 
ferred face centres and a strongly preferred phase of = tt. 
However, this is clearly a very special situation with a very 
particular resulting twist phase ip. Here we describe a more 
general test. The principle is to create a simple, asymmetric 
pattern in six faces of a spherical dodecahedron, translate 
them to the opposite faces, and apply a twist. By construc- 
tion, the cross-correlation should be very strong for this 
twist. In general, it should be weak for other twists. The 
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MCMC method is then used to try to detect the preferred 
twist. 

Some experimentation has indicated that the simplest 
test patterns are not necessarily sufficient to yield a con- 
vergent solution in 5-parameter space, since it is easy to 
induce false matches due to various periodicities between 
"adjacent" copies of the pattern. The function should be, 
in general, circularly asymmetric around each face centre 
so that only one favoured twist can be found. In the real 
Universe, it could happen that it is more difficult to find 
a single solution than in this artificial case. However, the 
result found in the WMAP data was that a single solution 
does indeed dominate. What is important to test is whether 
an artificial signal can lead to a strong, unique false detec- 
tion. 

The test pattern we adopted is 



W(100- 



(33) 



36000|6| ' 

where I and b are both in degrees. This is a smooth, mod- 
erately sinusoidal, asymmetric function of I and b, used for 
all pixels within 31.7° of the centre of a dodecahedral face, 
i.e. within a disc approximately touching the edges of the 
pentagonal border. These pixels are then copied by a trans- 
lation to the opposite face and rotated by a twist of 4>^. 

There is only a trivial algorithmic difference between 
stepping in all five parameters or stepping in a single pa- 
rameter. In this way, we can test most aspects of our 
method and software, as indicated in Sects 13.51 and 13.61 in 
particular, the estimates of the correlation functions, the 
calculation of the probability estimator defined in Eq. 
and the basic MCMC functionality, by starting our sim- 
ulated maps at the correct input values {l,b,0,a) for the 
simulated map and holding these constant, while starting 
at an arbitrary phase (f) unrelated to the input value of 
the simulated map. 

If something is wrong with the correlation function es- 
timation, the probability estimator defined in Eq. ((25)) or 
the basic MCMC functionality, or if a phase (p — ±7i'/5 
is somehow accidentally built into the software, then this 
should be revealed by toy simulations of this type. Figure[22l 
compares input and output values of (j) based on 20 MCMC 
chains, each of 6000 steps. Each input toy map is generated 
for a random point (I, 6, 9, a, 0), and the chain is started at 
the same b, 9, a) position but a random twist (j) £ [0, 27r]. 
Since the artificial maps contain matching 31.7° discs, we 
restrict the circle size in the MCMC chains to a < 30°. 

Estimates of 4>o are based on steps 1001 to 6000 (i.e. 
removing a burn-in of 1000 steps) of a simulation, using a 
minimum probability requirement of P > -Pjjiin = 0-5 as 
for the analysis on real data. The mean of the differences 
S = (j>Q — (p^ for the 20 chains is -1-8.3° with a standard error 
in the mean of 4.6°. Clearly, the input and output twists 
are statistically close to equal. Moreover, there is no sign 
that either of ±7r/5 are favoured due to either an intrinsic 
error in our method or to a software error. If the correla- 
tion function estimates, "probability" estimator, and gen- 
eral MCMC functionality contain any hidden errors, then 
these do not explain our estimate of ^ « +TT/5. 

5.4. Is the expected distribution of(j> uniform? 

If the simply connected model is assumed, then the spatial 
auto-correlation on large scales could conceivably be used 



to model the expected cross-correlations, though this would 
probably neither be easy nor free of model-dependent as- 
sumptions. 

For a fixed circle radius, as \4>\ increases from zero to tt, 
the separation between identified points increases. Thus, if 
the autocorrelation were monotonically decreasing as sep- 
aration increased, then higher values of \(j)\ would imply 
lower correlations. 

However, for circle radii a^ir/A, which constitute al- 
most all the points shown in Fig. [U the separation (as- 
suming simple connectedness) between opposite circles is a 
minimum of 2cos(7r/4) « 1.4 times an SLS radius, i.e. a 
minimum of 13.4/i~^ Gpc. If the auto-correlation function 

is reasonably isotropic, then over the range |(/)| G (0, tt), 
the cross-correlation will sample ^^((i >, 13.4/i~^ Gpc). 
Figure [T] shows that this function is approximately fiat 
and close to zero, at least relative to the correlations at 
d< 4/1-1 Gpc. 

The largest circle radii we considered here are those 
with a — 60°. This gives a minimum of one SLS radius 
(9.5/1^1 Gpc) separating a pair of points under the simply 
connected assumption, which is approximately the sepa- 
ration at which becomes flat and approximately zero. 
Hence, if the auto-correlation estimated under the simply 
connected assumption and shown in Fig.[l]is approximately 
correct, then the expected distribution of </> for the range 
of circle sizes considered here should be close to uniform. 

Independently of trying to explain our result in terms 
of a simply connected model, what other evidence could 
support the PDS interpretation? 

5.5. Can this solution be confirmed on smaller scales? 

Probably the biggest open question raised from the present 
work is the nature of the signal at smaller length/angular 
scales: r < 0.5 Gpc or 0^ < 3°. Investigating this would 
require much smaller Markov chain steps than we have used 
here. 

To have our Markov chains converge in a reasonably 
short time but also have them start at arbitrary positions 
in parameter space, we defined a 10° width Gaussian for 
stepping between points. Now that the globally optimal 
point in parameter space has been found, chains with a 
narrower width, and/or a stricter "probability" function 
for deciding how probable it is that a point in parameter 
space corresponds to a correct solution, could be used to 
find an optimal point with higher angular resolution, by 
starting from our present solution as an approximate solu- 
tion, rather than starting from arbitrary points in parame- 
ter space. 

However, if our solution is physically correct, then while 
there is no reason why the intrinsic density-density cross- 
correlation at sub-gigaparsec scales should differ from the 
density-density au^o-correlation, we can expect some differ- 
ences in temperature-temperature correlations due to vari- 
ous projection effects: both the Doppler effect and the inte- 
grated Sachs- Wolfe effect (ISW) should cause some differ- 
ences between the two functions when correlating tempera- 
ture fluctuations. Another complication on small scales, as 
Aurich et al. (2006) have shown, is that the WMAP ILC 
map has a lack of signal on these scales. 

Furthermore, could the stronger cross-correlation for 
our PDS solution at around 3-8h~^ Gpc relative to the 
auto-correlation be explained by ISW dampening of the 
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auto-correlation on these scales? Since the cross-correlation 
temperatures corresponding to these cross-correlations are 
observed, in general, at widely different angles, they are less 
likely to be foregrounded by the same ISW fluctuation, so 
this seems possible at least qualitatively. 

A practical problem in extending the present method 
to smaller scales is the need to use enough pairs of points 
(xi-^^Xi^) on the SLS for which d{xi-^, [9ji^i2)]) is "small". 
For a given, finite number of points TVp, there are fewer 
and fewer pairs of points separated by smaller and smaller 
separations. 

As explained in Sect. 13.51 we used iVp = 2000 points dis- 
tributed uniformly on the sphere. In fact, a larger number of 
points are generated when a galactic contamination mask is 
used. To avoid any chance of bias, the points are generated 
uniformly, excluding those falling within the mask, until the 
required number of points have been obtained. Excluding 
distantly separated pairs of points, where the distance is the 
minimum of {^(a;^^, [gjixi^)])} for the 12 holonomy trans- 
formations gj, is a heavy computational task equivalent to 
making the main calculation loop needed for correlation 
function estimation. 

Testing of smaller scales therefore requires considerably 
more pairs of points, and this would seem to scale as iVp, in- 
creasing the calculation time proportionally. Improvements 
to this scaling using geometrical arguments would likely 
be complex and would risk inducing an ambiguity as to 
whether the new results are due to the non-uniform dis- 
tribution of the points or to the intrinsic properties of the 
data. 

5.6. Are the "secondary" and "tertiary" solutions explainable 
as PDS harmonics (eigenmodes)? 

If the eigenmodes of the PDS imply features such as those 
shown in Fig. [5] and Table [31 then this would support the 
PDS interpretation of the signal found here. 

5. 7. Is this solution compatible with other PDS analyses of 
the WMAP data? 

The Roukema et al. (2004) "hint" of a solution was based on 
looking for identified circles, so, as noted above in Eq. (O, 
was based on using a much smaller part of the available 
signal to search for a candidate favoured PDS orientation. 
Key et al. (2007) and Lew & Roukema (2008) agree that 
the solution is a local maximum in the parameter space, 
but that it is not statistically significant. Using the greater 
amount of information available from the cross-correlation 
at separations going up to r < 4/i~^ Gpc rather than just 
r = 0, the results presented in Sect. 14.41 agree that the 
Roukema et al. (2004) PDS solution is not globally favoured 
by the data. Moreover, MCMC chains starting at that solu- 
tion move towards and converge towards the solution found 
with the present method. 

One way of comparing with other PDS analyses, which 
use quite different methodologies to analyse the WMAP 
data, is to consider estimates of the total density param- 
eter r^tot- Luminet et al. (2003) initially favoured an esti- 
mate of i^tot ~ 1.013, but in more recent work using many 
more eigenmodes of the PDS, they favour iltot ~ 1.018 
(for rim ~ 0.27) by requiring maximal suppression of the 
quadrupole (Caillerie et al. 2007). In analyses based on 



WMAP first-year C/ estimates, Aurich et al. (2005b) found 
1.016 < Otot < 1.020, and using a combined match circles 
method, Aurich et al. (2006) favoured iltot ~ 1.015 (for 
nm « 0.28). 

Are these consistent with our solution? The total den- 
sity parameter f^tot relates tightly to the circle radius a, 
with some dependence on ilm, e.g. Fig. 2 of Roukema et al. 
(2004) or Fig. 10 of Aurich et al. (2005b). The estimates 
n^ot = 1-015 (for = 0.28) and ^tot = 1-018 (for 
rim = 0.27) correspond to a = 40.5° and a = 47.6° respec- 
tively. Gundermann (2005) also found il^ot — 1.017 (for 
= 0.26) as one of his two solutions; this corresponds to 
a = 46.8°. 

As discussed in Sect. 14. 2[ the "matched discs" inter- 
pretation (Fig. [S]) of Fig. [5] would imply a true matched 
circle radius of ~30-40°. This is clearly compatible with 
the Aurich et al. (2006) estimate, and possibly compatible 
with the Caillerie et al. (2007) and Gundermann (2005) esti- 
mates. Clearly, a variety of different methods, either requir- 
ing suppression of large-scale power (other work) or max- 
imising large angular scale correlations in a PDS-like sym- 
metry (our work), and either with the in-built constraint 
that the twist is ±7r/5 (other work) or without that con- 
straint (our work), seem to converge on similar results. 

6. Conclusion 

It seems hard to avoid the conclusion that cross-correlations 
of temperature fiuctuations on would-be adjacent copies of 
the SLS, which are distant from one another and are (on 
average) very weakly correlated according to the WMAP 3- 
year observations, (i) imply a highly cross-correlated "gen- 
eralised" Poincare dodecahedral space symmetry (Figs [S]- 
[3 Table [1]) at which these on-average uncorrelated fiuc- 
tuations happen to be well-correlated with one another 
fFig. [TT|). and, moreover, (ii) this favoured solution is dom- 
inated by a signal whose twist phase (j) lies within a few 
degrees of one of the two twist phases necessary for a valid 
PDS model (Tabled. 

These two successful predictions of the PDS model fol- 
low the WMAP confirmation of the generic prediction of 
small universe models, a power cutoff in structure statistics 
on large scales, and the solution appears to be consistent 
with quite different PDS analyses of the WMAP data. 

Do we really live in a Poincare dodecahedral space? 
Further constraints either for or against the model are cer- 
tainly still needed, but the evidence in favour of a PDS-like 
signal in the WMAP data does seem to be accumulating. 
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